Asteroids, Comets, Meteors 1991, pp. 359-362 
Lunar and Planetary Institute, Houston, 1992 


359 

Numerical Simulations of Cometary Dust j /y p ^ 4/ 0 

N 9 3^“19 1 9 3 


David J. Lien 
Bucknell University 


Most observations of comets are done photometrically or spectrophotometrically. The 
interpretation of the aperature-averaged flux is relatively simple for an isotropic, radially 
expanding coma of infinite extent - the canonical model. However, the interpretation of 
the observations is not so clear when the the motion of the dust is affected by radiation 
pressure, or when the emission is time-varying and anisotropic. For example, in a sample 
of CCD images of 10 comets, Jewitt and Meech (1987, Ap.J. 317, 992) found that the 
photometric profiles of only three comets were consistent, within the observational errors, 
with the profiles predicted from the canonical model. Photometric observations with large 
apertures, however, seem to suggest that the canonical model may be quite adequate (c.f. 
Osip, Schleicher, and Millis, 1992, Icarus 98 , 115). 

The dust itself is characterized by a size distribution, with size dependencies on the 
expansion velocity, the scattered and thermal radiation, the response to radiation pres- 
sure, and probably the density. How good then axe the approximations normally used in 
determining the production rates of the dust when these effects axe present? 

Given an isotropically expanding dust envelope (from a non-moving nucleus), the 
radial profile can be shown to decrease as the inverse of the projected cometcentric distance: 
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where a(p ) is the column density at the radial distance p, Q is the production rate [s 1 ], 
and v is the expansion velocity. If the limits are not ±oo, then 
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where R is the edge of the envelope. 

The total number of particles inside an aperture of projected radius s is 
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where the last equality is for R = oo. Since the time over which the dust is emitted is 
finite, the total number of dust particles in the coma is just the product of the production 
rate and the interval of time over which the dust was emitted. Evaluating equation 3 for 
an aperture equal to the outer edge of the coma yields N(R) = = Qt , as expected. 

As paxt of a program to better understand the dynamics of cometary dust and gas, a 
computer program has been developed which numerically simulates the emission of both 
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dust and gas from the nucleus of a tilted rotating nucleus. Only the dust coma will be 
discussed here. The dust coma simulation includes the effects of the size dependencies on 
the expansion velocity, and scattering or thermal emissivity (based on either approxima- 
tions or Mie theory calculations) and on the response to radiation pressure. Anisotropic 
emission is approximated by a gaussian jet centered at any latitude and longitude on a 
rotating nucleus of arbitrary rotation rate and obliquity. 

The “image” of the gas or dust coma can be generated, as well as aperture- or annulus- 
averaged fluxes. An example of the annulus-averaged flux is presented below for P/Halley 
on 15 March 1986. Figure la shows a(p) from the numerical simulation of dust emission 
for 14 days at v — 0.025 km/s with j3 = 0. This gives a maximum size of the coma of 
about 44”. A least squares fit of equation (2) to the data is displayed as the solid line in 
Figure la. The solid line in this and subsequent figures corresponds to the theoretical fit, 
and the dotted line is derived from the numerical simulations. The aperture averaged flux 
is shown in Figure lb, along with a least squares fit of equation (3) to the data. 

Since the production rate is constant and known from the parameters input into the 
numerical simulation, the R = co solution can be solved for Q based on the results in 
Figure lb. This is shown graphically in Figure lc, where the production rate, Q , derived 
by assuming the R = oo solution to equation (3) is plotted as a function of the aperture 
size, s. The thin horizontal line indicates the true production rate. For this simple model, 
the R = oo approximation overestimates the production rate for small apertures, and 
underestimates Q for large apertures. As expected, the agreement between the R = oo 
approximation and the simulations inproves as R increases. 

The dynamics of a large dust particle once well away from the nucleus depend primarily 
on the initial velocity vector of the particle and the ratio of the radiative pressure force to 
the force of gravity, (3. Small charged dust grains are not considered here. As /? increases, 
the coma becomes increasingly non-spherical, and one expects that the application of 
spherical approximations to the coma will become increasingly inaccurate. To test this 
statement, a series of models were run with v = 0.025 km/s and a single (3 from the 
range 0.0001 < (3 < .5. These parameters were chosen such that a spherically expanding, 
(3 = 0 envelope would occupy approximately ^ the length of the image array. Figure 2 
(/ 3 = 0.0001) and Figure 3 (/3 = 0.1) show both a projection against the sky of a subset of 
the total number of particles used in each simulation (typically 50,000), and the value of Q 
as a function of the aperture diameter. Again, the thin horizontal line indicates the input 
value of Q In these models, even a small (3 causes a significant change, with the predicted Q 
more than a factor of 2 too low for most aperture sizes. It is interesting to note that these 
simulations suggest that spectophotometric observations, with relatively small apertures, 
will yield higher dust production rates than the larger aperture photometric observations. 

The question of a size distribution on the predicted Q is illustrated in Figure 4. 
To illustrate the effects of a size distribution, a (3 distribution which is rectangular in 
log((3 ) (equal numbers of particles in equal size bins in log(/?) with /3 within the range 
0.0001 < /3 < 0.1 was used. Each particle contributes an equal amount to the flux, and 
the velocity scales as yfft. The scale factor is chosen so that the velocity of the center of 
the (3 distribution is 0.025 km/s. Figure 4 shows the results of this simulation. 
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To more realistically simulate the scattering or thermal emission properties of real 
grains, the same input parameters used in Figure 4 were modified by having each parti- 
cle contribute a flux proportional to (3 ~ 2 . This relationship is derived from the inverse 
relationship between the diameter of the particle and the value for /?, which is valid for 
particles large compared with the wavelength of light. The results are shown in Figure 5. 
To determine the production rate, an “average” flux and velocity must be chosen. The 
values for the center of the /? distribution was used to determine Q from the simulations. 
The production rate is overestimated by as much as 2 orders of magnitude ( Q = 0.041 s -1 ), 
even with a large aperture. This result suggests that the derived production rate is very 
sensitive to the the choice of the optical properties of the “average” particle. 

Although the results of most of the simulations run with a single beta are consistent 
with the statement that the production rate derived by assuming R = oo is correct to better 
than a factor of 5, this conclusion is misleading. As illustrated in Figure 5, when a broad 
size distribution with the intensity proportional to f3 ~ 2 is used, the derived production 
rate is off by orders of magnitude. The main difficulty lies in determining the “average” 
velocity and /3. Future work will concentrate on more realistic models: longer emission 
time, incorporation of the optical properties of spheres composed of potential cometary- 
dust materials, and a more realistic size or /? distribution. Additionally, efforts will be 
focused on using the results of these simulations to aid in the interpretation of photometric 
and spectrophotometric observations. The results from the simulations presented here 
suggest that the canonical model, which assumes a radial outflow of dust, is inadequate 
in measuring the dust production rate. However, further simulations with more realistic 
input parameters are needed before any final conclusions can be drawn. 
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Figure 2 





